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Abstract 

The Landau expansion for the free energy of the superconducting mixed state near the upper critical 
field in powers of the square modulus of the order parameter averaged over Abrikosov lattice is derived. The 
analytical calculations has been carried out in frame of Gor'kov formalism for 3-dimensional isotropic BCS 
model beyond the limits of quasiclassical approximation, another words with Landau quantization of the 
quasiparticle energy levels taken into consideration. The derivation is performed at low temperature and 
high enough but finite crystal's purity. The effect of Pauli paramagnetic terms is taken into account. The 
quantum oscillations of the critical temperature, the order parameter's amplitude and the magnetization (de 
Haas-van Alphen effect) in the mixed state are found. The limitation of validity of a mean field approach 
due to critical fluctuations (Ginzburg criterion) for the phase transition under consideration is established. 

1 Introduction 

The growing body of measurements (see the most recent publications HJ, l^j,^, §|) demonstrates the de Haas - van 
Alphen effect in the mixed state of type-II superconductors. The theoretical description of this phenomenon is 
rather cumbersome and as a result we have a vast amount of publications devoted to the subject [||, ^, 0, ||, ||, 

origination of theoretical troubles is a nondiagonality of the order parameter matrix 



in Landau representation preventing a correct derivation of the quasiparticle energy spectrum in the mixed 



state. The later can be calculated properly just near H& [13]. The intention to avoid of this problem leads to 
idea to develope a theory far below H C 2 where a vortex core radius is much smaller than the distance between 
vortices and one can assume the space constancy of the order parameter modulus. This region on the first sight 
intersects with field interval of the observations of the dHvA oscillations in the mixed state extending down to 
the fields of the order of 1/2 or even 1/5 of the upper critical field (see Q). However in the superconductors 
with very large Ginzburg-Landau parameter k, which is of the order of 20-30 in typical for observation of dHvA 
in the mixed state materials, the field interval of observation should still be considered as relating to vicinity 
of the upper critical field where the distance between the vortices almost coincides with the superconducting 
coherence length or vortex core radius. The former starts to exceed the later only at the fields of the order of 
H&/k. So, in fact we need in a theoretical description of dHvA effect in a vicinity of the upper critical field 
(H C 2 — H) /H C 2 < 1 where one can assume the field inside of superconductor as coinciding with homogeneous 
external field and use the Abrikosov ansatz for the periodic in space order parameter distribution. 

Here we present the derivation of the Landau expansion of the free energy density in vicinity of the upper 
critical field in powers of the square modulus of the order parameter averaged over Abrikosov lattice, 

F s -F n = aA 2 + ^A\ (1) 

which is valid until (H C 2 — H)/H C 2 < 1 in spite of relative value of the order parameter amplitude and the 
distance between Landau levels (for more exact condition of validity (|l|) see (p7|) ). The treatment of the problem 
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is carried out in frame of the Gor'kov formalism with Landau quantization taken into consideration. It has 



been done earlier for 2-dimensional electron gas either numerically ||, 10, 11] or analytically but out of dHvA 
observation region |l2|] when the temperature T exceeds the distance between the Landau levels [] uj c = H/m*. 
The reason for the treatment of the problem in the unsuitable region is divergencies of a and (3 coefficients 
in the formula (|lj) at temperature T — > in absolutely clean limit. Unlike to these papers we shall treat the 
problem analytically for 3-dimensional isotropic BCS model and escape of the divergencies by the introducing 
a width of Landau level V originating of impurity scattering or a specimen's inhomogeneity. In addition to 
1 13] where a and (5 coefficients has been determined neglecting the oscillating terms here we take them into 
account. The expressions for the functions describing the oscillating behaviour of the critical temperature, the 
order parameter's amplitude and the magnetization in the mixed state are found. As in a normal metal the 
tiny oscillating corrections to the free energy are transformed into the big summands in the magnetic moment 
as the result of differentiation of very fast oscillating functions of the ratio of Fermi energy to the cyclotron 
frequency. It is shown that at the distance of the order of £f/T c oscillations from the upper critical field there 
is noticable decrease (of the order of unity) of the amplitude of dHvA oscillations of the magnetization in the 
mixed state. 

As was pointed out in the paper [14] the Landau level splitting due to Pauli paramagnetic interaction of 
electron spin with external magnetic field can cause the change of character of the phase transition to the 
mixed state. Namely, in clean enough type-II superconductor the second order transition transforms into the 
first order transition at temperatures below some tricritical point depending of impurity concentration and 
determined by equation @(H, T) = 0. Here we confirm that there is a principal opportunity of such a kind 
phase diagramme at low enough temperatures T < T when the Zeeman splitting of Landau levels exceeds 
a width of them: [i e H > T. However in view of very complicated analytical structure of the expression for 
f3(H, T) it is difficult to determine definetely its sign in the region mentioned above. One may assert only that 
at fj, e H < r the second order type transition takes place up to zero temperature. 

The investigation of the fluctuation effects is also undertaken in the present work. Unlike to the zero field 
case where the fluctuation region in vicinity of the critical temperature is of the order of {T c /ef) a that is 
negligibly small, the critical region at low temperatures in vicinity of the upper critical field is of the order 
of (Tc/ef) 2 . The period of critical temperature oscillations has the same order of magnitude. This means 
that the field interval for the reentrant normal metal-superconductor phase transitions originating of critical 
temperature oscillations is located in the critical region. Hence the reentrant behavior has the poor chances 
for observation. It is shown however that the quantum oscillations of T C (H) can be accessible for observation 
in lowest temperature region T < T <C u> c . 

The paper is organized as follows. The formula (|l]) is derived in the next Section. The oscillating with 
magnetic field behavior of the a (or another words the oscillating behavior of the critical temperature) and j3 
coefficients are described in the following two Sections. The magnetizations oscillations in the superconducting 
mixed state are calculated in Section 5. The critical fluctuations are discussed in the Section 6. 



2 Free energy 

The difference of the superconductor - normal metal free energy densities expanded in powers of the order 
parameter A(r) and averaged over a superconductor volume V has the form 

F s -F n = ^J dr(F s (r) - F n ) = J dr|A(r)| 2 -yj dr 1 dv 2 K 2 (r 1 ,r 2 )A*(r 1 ))A(r 2 )+ 

— y , dridr 2 r3dr 4 K4(ri,r 2 ,r3,r4)A*(ri)A(r 2 )A*(r3)A(r4). (2) 

Here g is the constant of the pairing interaction. Neglecting all the vortex corrections that is valid in ultra 
clean limit T < co c considered here we shall use for the functions K 2 and K± the following expressions 

K 2 { Yi ,t 2 ) = \ T^G CT ( ri ,r 2 ,^)G- CT ( ri ,r 2 ,-^), (3) 

1 o=±\ v 



1 Planck's constant h, electron's charge |e| and velocity of light c are chosen to be unity throughout the paper. 
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^ 4 ( ri ,r 2 .r 3 ,r 4 ) = \ £ G CT ( ri , r 2j a^G^fo, r 2 , -^)G CT (r 3 , r 4 , c^G^n, r 4 , -w„), (4) 
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Z (7=±1 f 

ccv = uj v + Tsign uj u , u> u = irT(2v + 1), 

written by means of the electron Green functions G a (ri, r 2 , lo u ) in normal state under an external magnetic 
field. The later are related with the Green functions depending only on the relative coordinate G a (ri — r 2 ,u) l /) 
as 

G a (ri,T2,w v ) = exp U J A(s)dsj G a (ri - r 2 ,cD„). (5) 
Making use this formula one can rewright the quadratic in A(r) terms in the formula @ as 

(I- J dRexp [~ l -He 2 ^ tf 2 (R)) \ J dr\A(r)\ 2 , (6) 

where g 2 = R 2 — Z 2 and 

K 2 (R) = \Ys T^G ct (R,^)G- ct (R,-^)). (7) 

a=±l v 

As it was discussed in Introduction for the strong type-II superconducting materials the Abrikosov solution 
for the order parameter 

A(r) = A/(r), (8) 

v=integer \ J 

is valid in broad enough vicinity of the upper critical field. Here for simplicity we have chosen the Abrikosov 
square lattice solution in the Landau gauge A(r) = (0, Hx, 0). The elementary cell edge length a is such that 
a 2 = tt\ 2 , where A = H^ 1 ! 2 is magnetic length. We assume that magnetic field is uniform and coincides with 
an external field. This is certainly true in vicinity of H C 2 for the superconductors with large Ginzburg-Landau 
parameter . 

Using formulas (||) one can rewrig ht (|) in the shape ([[]) 

F s -F n = aA 2 + ^A\ 



where 

a 



--J dKexp (-7^£ 2 ) K 2 (R) (10) 



g 

and 

&=\\ *idr 2 r 3 dr4K4(ri,r 2 ,r3,r4)r(r 1 )/(r 2 )/*(r 3 )/(r 4 ). (11) 

The amplitude of the order parameter A has to be found from the free energy minimum condition. 

The a and (3 coefficients are the oscillating functions of the magnetic field which will be obtained in the 
following two Sections. 

3 Critical temperature oscillations 

The critical temperature as a function of magnetic field at T <C T c have been found in the paper by L.Gruenberg 
and L.Guenther |[(|. It is determined by the equation 

a(H,T)=0, (12) 

where 

Nn { H — H c2o T 2 3/2 (lO c \ 1 / 2 o/n I/O W c 



T ^ 1 _ L ^ L + J_ _ ^S/2 p Si _ 23/2^1/2^^ (13) 

2 H c2o T£ V H J M I 
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ft = £ f ex P p™(^ + ^/m (14) 



n=l ^=0 



Si = —K±± E(-D m+ "exp f-^ + n ^ + ^) C ° S ( 2 ; |n " m 'f"^, (15) 
^ „=i™>„m V ^ c 7 |n-m|V2 

H c2o = ^--f . (16) 

Here, /x is chemical potential, t>i? is Fermi velocity, T is Landau level width, [i e is the electron's magnetic 
moment, Nq = m * kp/2ir 2 is the normal metal electron density of states per one spin projection. S*o is slow 
(logariphm) function of magnetic field with numerical value of the order of unity. 

The substitution of the first two terms of (|l3|) in the equation ([H|) reproduces of the quasiclassical result 
found by L.Gor'kov |7] 

H = H c2o (l - ^S ) (17) 

for the upper critical field at temperature tending to zero . The third and the fourth terms in the formula ( Jl3 ) 
give correspondingly nonoscillating and oscillating corrections to Gor'kov's expression owing to quantization of 
the quasiparticle levels. For the further considerations it will be convinient to represent the functions S\ and 
S2 in more simple analytical form. 

The summation over v in (14) yields 

- T 00 expf-^M) 

Si = —&y \ 2 V , (is) 

UJ C ^ sinh^^ ^ ' 

II — 1 (jj c 

where 

a = T + ir)u c , (19) 

and 

i]lo c = /j, e H - -oj c , (20) 

is the Zeeman energy deviation from a half integer number q/2 of distances uj c between Landau levels. For 
noninteracting electron gas in a perfect crystal a = and 

*"S ta iSr (21) 

diverges at T — > 0. This result formally means an existance of superconductivity of noninteracting electron gas 
in a perfect crystal at T = in arbitrary large field |l(| . 
In the opposite limit ttT <C \a\ when 



n=l 

At \a\ <C co c /4it it is 



47T 47r|a| 24 Via 
For 4-7rr ~ oj c one can estimate sum (O) by its first term 



Sx = Un^h-^[h) • (23) 



ttT ( A.v;T\ 

Sl = . ; 47T 2 T exp cos Attt]. (24) 

u r sinn L V w c / 



The substitution ( |24| ) into equations (|13|), ( |12| ) provides us by the averaged over oscillations behavior of upper 
critical field at low temperatures 

H c2 {T) = Hc Jl + 1 -^(l- 8 4(l-) 2 ) exp cos 4^1 (25) 



/i \ 3 \Ldc2J j V Wc2 
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Here uj C 2 = H C 2/m* is the cyclotron frequency at the upper critical field. Let us denote the number of Landau 
levels at H = disposed below the Fermi level ef as 



uj c 2tt 2 V T c 



The number n C 2 in typical materials with small Fermi energy and relatively high T c where dHvA effect in the 
mixed state has been observed is of the order of one hundred. So, the upper critical field at zero temperature 
is noticably larger (~ n c2 ) than its quasiclassical value (|i~6|). The decreasing H C 2 with temperature follows 
T 2 law. This parabolic dependence in the region 

2vr 2 E F 

is much faster than Gorkov's parabola (|l7|). 

For 5*2 the summation over v in (14) produces 

t °° 00 exo( Mm+n)(r+i^ e H) \ cqs ( 27T ( m n \j± A 

Changing the summation variables to n and m — n = I we get 

vrT ~ (-1)' / MW)\ / n vr^ exp(- ^/ ^ 

52 = — exp I — ^ — ) cos - 4 J £ sinh ^(^ • (28) 

Similar to S\ for noninteracting electron gas in a perfect crystal ((2n e H = uj c , T = 0)) the expression for S2 is 
diverged at T — ► (see also [|16| ) . 

In more realistic case 47rr ~ w c one can keep only the first term in the sum over n: 

7r r^ (-l)'exp(-^g± 2 l) / ^ /^l[(J + 2)\ r9Q ^ 

Do = > » on , oW COS Z7T/ COS . (29) 



Substitution of (24) and (^) into ( |T2| ) and the resolution of the equation a(//, T) = in respect of T gives 
the oscillating behavior of critical temperature. The simple estimation shows that the amplitude of oscillations 
at T ~ u> C 2 ~ 27r 2 T c 2 fep is of the order of 

^ ~ (I)' (30) 

At lower temperatures the amplitude of oscillations is somethat larger. The period of one oscillation 8H is easy 
to find from 

I - $ n - §EL - £p 6H 
uj c uj c H 

So, 

SH o fT c \ 2 

As we shall see later this value coincides with the critical region in vicinity of upper critical field at low 
temperatures. Hence the critical temperature oscillating behavior as a function of magnetic field exists only 
in frames of mean field approximation and has poor chances to be observable. On the other hand outside the 
narrow region of strong critical fluctuations in the mixed superconducting state the function 

a{H,T) = a(H,T) + a osc (H,T) (32) 

consists of smooth function 

<x{H,T) = — (33) 

^ tl c 2o 
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(see (p5h) and fast oscillating function of magnetic field 



a osc 



,rr T , i/2 3/2 6XP (~ gEl £ ±gl ) f, ,P *\ { 2TT» e H(l + 2) \ 
(H,T) = -2' ixl N -^ - ^ cos ( 2ttZ— - - cos (34) 



V/sinh^±^ V w e 4 

giving the main contribution to the oscillating part of the magnetization. 

4 Oscillations of (3(H, T) 

Instead of formula (pi]) it is more suitable to use an expression for in the so called magnetic sublattices 
representation (see |I~3|]) 



(3{H,T)= l -Y. T Y. E \^G-°{i nk ^-u v )G°(£ mkz ^ v )G-°^ (35) 

<r=±l v n,n' ,m,m' 

where 

G°(H nkz ,^) = ~rz j^— = , (36) 

1 k 2 

ink z = u c {n + -) + 7-^ - M> (37) 

2 2m* 

-^mm' = y ■^^2/nm(q)/^ m '(q)/n'm(q)/n'm'(q) (38) 

and / nm (q) are the matrix elements of functions (|9|). 

In view of extremely complicated structure of general formula (|35| ) we are limited ourselves by the analysis 
of diagonal terms n = n' = m = m' giving the main contribution to the expression (|35"|) at T < u c . For this 



case as was shown in the paper |L3| F™ rs 1/(2tt\) n and 



/3(jff ' r) = 2f2^AW ^ ^(G°(U Z ,^)G-°(U Z ,-^)) 2 . (39) 

' ct=±1 f n=0 

Let us represent 

(3(H, T) = p(H, T) + f3 osc (H, T) (40) 
as the sum of smooth nonoscillating function and fast oscillating function of magnetic field. For the calculation 



of nonoscillating part of f3 following the papers [L4[ one can substitute the summation over n by the 
integration according to 

1 ~ r dk 7 r°° dn f dk 7 „ r _ W 2 



2ttA 2 n 

n=0 

Making use n ~ //sin 2 #/u; c and performing the integration over 9 we get 

- = £ r^iv / di (g*(z,l,)g-°&-u,)) 2 , (4i) 

The integration over £ yields 

mr) . Wo ^ T ££zM!, (42) 

This expression can be positive or negative depending on the relation between the values of T and [i e K . At 
T = K when it is possible to change the summation to the integration we obtain 

- „ "cln£ r 2 -Qu e ff) 2 

^ (jff ' 0) - ^ 16^ c^ + C/i^W ( 43 ) 
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This result demonstrates the positiveness of (3 up to zero temperature at V > fJ, e H. In the opposite case 
(3 at T = K is negative. Hence below some temperature ~ [i e H the second order type transition to the 
superconducting state starts to be of the first order. Becouse we have thrown out the nondiagonal terms it is 
difficult to say definetely is it the case or not and if it is so what is the temperature of tricritical point. Leaving 
this problem for future experimental and theoretical investigations let us assume positive value of 

which is valid at least at the fulfillment of the conditions T, fi e H < T < u c . 

To find the oscillating part of P(H, T) at the same conditions let us apply to the expression ( |37| ) the Poisson 
summation formula 

[3 osc (H, T) = 

here e = Cn,k z + M- Subsequent calculations follows by well known Lifshits, Kosevich derivation [18] of normal 
metal de Haas - van Alphen oscillations. Performing the saddle point integration over k z we get 



exp [2ml 



d£,jrz — m v 2 , .- c/ c ir\2' ( 46 ) 

The integration over £ yields 

0°° C (H, T) = N 0J ^ Ez=i cos ( 2 vrZ£ - f ) «T f) exp (-2 ^ + ^ ) 
At the fulfilment of the conditions T, fi e H <T < uj c one can rewrite this expression in more simple form: 



3/2 



where 



/( r) " Z (y+l)3 + (y+!)2 ^'^I> (-■'•//)• 







Thus, as it is expected, the oscillating part of (3 (|47|), ( |4q ) is \/^c2 times smaller than its smoth nonoscillating 



part (H), |_ 

The preliminary calculations of /3 has been published in IflQf] . Although the results of that paper are in 
general qualitatively correct there there are the errors in the values of j3 and (3 osc (compare formula (^) in [|0|] 
and (fl||), (fi8|) in the present article). 

5 Magnetizations oscillations 

Minimization of the expression (|lj) over A yields 
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Now we can calculate the magnetization 



M, 



dF s 
dH 



Keeping at the differentiation only the fast oscillating terms we get 



ada osc a 2 dp osc 



(3 dH 2~p 2 dH ' 



Here M° sc is the oscillating part of normal metal magnetization [18| including Dingle factor exp(- 
due to scattering by imperfections 



M osc 



1=1 



Here 



2V 2 nH \uj Z 3 / 2 sinhA/ V ^ c 4/ V tu c J \ w c 



A, 



2^ 2 /T 



For the subsequent terms in ( |5l| ) making use (|33|), ( |34|) and approximative expressions (Q) and 
smooth and oscillating parts of P we get 



(50) 

(51) 
-2nW/uj c ) 

(52) 

(53) 

(54) 
for the 



a da 



p dH 



27/2^2 ^ Hc2 (T)-H 

/i 7T ' 



E h m o\ ■ h x sin 2nl ~ I cos 2?r z + 2 exp - 2n z + 2 - 

^ (/ + 2) smh Az+2 V w c 4 / V w c / \ w c 



(55) 



a 2 8f3 c 



N - 



2 5/2 7r2 p2 



3/2 



'H c2 (T)-H 



2j3 2 dH H 

OO / 

^(_l)'+l/l/2 sin ( 2 7T/ 



3/2 In 2 JL 



1=1 



H c 2o 

r 



(56) 



The corrections to the dHvA amplitude in the mixed superconducting state have the same structure as the 
normal metal part: each term in the sums (|55|), (|56|) oscillates with the frequency lfi/u> c corresponding to the 
extremal crossection of the Fermi surface. From these expressions at 47rr ~ uj c it is easy to see that until 



H c2 (T) - H 



H 



< 



c2o 



the following inequality takes place 



\Mi 



> 



P dH 



> 



— in — 



a 2 dp° 



2P 2 dH 



(57) 



(58) 



So, in the limit of validity of the condition (1571) one can neglect of the third term in the formula (JsTT) and 



treat the (55) as the main mixed state correction to the oscillating part of magnetization in the normal state. 
Summing up (15^) and (|55|) we obtain 



rv 8a osc J5. 

M osc = M osc + = £ MlMgh 



1=1 



M sl = l-^ 



l_/4vrr\ 2 fi 1 / 2 H c2 (T)-H Z 2 A i+2 sinhA z 

7T \ U c J 



COS 



(2tt(Z + 2) 



1/2 , U 



H c2o (I + 2) A/ sinh A /+2 cos ( 2t\I^L 



47rr\ 
( ' X T ( — —) 



(59) 



(60) 
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This expression presents the main result of the paper. Looking on ( |60| ) one can conclude that at fi e H ~ co c /2 
as well as at fi e H <C uj c /2 the amplitude of de Haas - van Alphen effect in the superconducting mixed state is 
less than in the normal state. For the intermediate values of Zeeman splitting like fj, e H ~ w c /4 the amplitude 
of the magnetization oscillations in the mixed state can be even larger than in the normal state. 

The region of validity of the result ( |59| ) is determined by the inequality (57) which can be rewritten as (see 

©) 

H c2 (T)-H < Inn* (gi) 



H C 2o \fna. 
If we remember that period of oscillations (|3l|) is of the order of 



1 (62) 



H c 2o n C 2 

that means the formula ( |59| ) describes the oscillating part of magnetization in the field interval corresponding 
to 

1 /n^2lnn c2 . (63) 

oscillations below H C 2- This number in the typical for observation dHvA effect in the superconducting mixed 
state materials is of the order of several tens. Out this field interval the corrections to the normal metal 
amplitude of any order ~ (H C 2 — H) n shall be the same order of magnitude and the result (^) is inapplicapable 
to the mixed state dHvA effect description. It is worth to be noted here that the region (|63|) coincides with 



the region of the existance of gapless superconductivity found in the paper p3 |. So, out the field interval ( |571 ) 
below upper critical field the de Haas van Alphen effect in the superconducting mixed state has to be negligibly 
small. 



6 Critical fluctuations 

The problem of critical behavior of type II superconductors near the upper critical field had been considered 
in many theoretical papers (see for instance [^] and references therein). We shall not be interested here by 
the solution of this problem as whole but just an estimation of a width of the fluctuational region at very 
low temperatures in very clean materials. To establish the region of importance of critical fluctuations under 
magnetic field near the zero temperature let us remind first how to estimate it in zero magnetic field near 
T c . The most convinient for our purposes way is to compare the density of energy of critical fluctuations 
Ffi ~ T/(£(T)) 3 with mean field energy density F m f ~ Nq(T — T c ) 2 above critical temperature. Here 



t(T) = Co 



T-T c 



is the coherence length. The comparison of these values leads us to well known condition (Ginzburg criterion) 

^ < m 

of importance of critical fluctuations. To get the similar condition near H c2 at low temperatures one needs 
to take into account the one dimentional character of fluctuations along magnetic field direction such that the 
density of energy of critical fluctuations will be 

Ffl ~ 2^{H) ' (65) 

Where A = H^ 1 ^ 2 is magnetic length coinciding at H = H c2 (T = 0) with zero temperature coherence length 
Co = v f/2ttT c . The field dependent coherence length above H c2 is given by 
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Here we have omitted the logariphmic corrections to £(H) (see 
The density of mean field free energy has the form 

« 2 r, Ar £ fT 2 (H-H c 



The comparison of equations ( p5| ) and ( p7\ ) shows that the critical fluctuations contribution to the thermo- 
dynamical values surpass the corresponding mean field values at 

This inequality establishes the region of importance of critical fluctuations (Ginzburg criterion) above H C 2 at 
low temperatures. The applications of this criterion has its own limitations. The point is that it is valid 
for classical or thermal fluctuations with energies ~ T larger than minimal frequency of the order parameter 
fluctuations 

—fT^ < T ■ (69) 
ti c 2 J-c 

When temperature decreases, the region of importance of thermal flugtuations ( |68| ) shrinks slower (oc T 2 / 3 ) 
than the region ( |69| ) of applicability of Ginzburg criterion. So, at lowest temperatures the formula ( |68| ) is 
inapplicable. Let us compare vicinities fl68| ) and ( |69| ) say at 

UJ c2 7rT c 2 



/i 



4vr 2e F 

For the region of importance of critical fluctuations we get 



(70) 



< UvrMn^ ^ (71) 

H c2 \ UJ C J \£ F J 



and the thermal fluctuations region (|69| ) is limited by 

H — H c2 it T c 

—n-^ < (72) 

H C 2 I £f 

We see that the region (71) is still parametrically more narrow than the region (72). However due to the large 
numerical factor in (|7l| ) these regions can be proved of the same order of magnitude. 

At lower temperatures we have to calculate the quantum fluctuation contribution. For the problem under 
discussion this contribution is exponentially small, so one can say that the region of the lowest temperatures 
is fluctuationless. 

Turning back to the inequality (^) and taking into account the estimation (0) one can conclude that the 
region of importance of thermal fluctuations at least for temperatures above T ~ T presented by ( |70[ ) is given 
by inequality ( p8| ) that is broader than the period of critical temperature oscillations (|3l|): 

SH _, n _ 2 (T c \ 2 



^ — 

H c2 \e F ) 

Thus in this temperature region the oscillations of critical temperature have poor chances for observations. On 
the other hand at temperatures 

T < T < ^ (73) 
we have region free of fluctuations where one can hope to find the quantum oscillations of critical temperature. 



7 Conclusion 

The quantitative theory presented in the paper demonstrates the suppression of the de Haas van Alphen effect 
in the superconducting mixed state according the equations (|59|), (|60[). The suppression is developed in the 
region below the upper critical field embracing up to hundred oscillations of magnetization ( |63| ) and limited by 
the inequality (|57[). The quantum oscillatons of the critical temperature are accessible for observation at very 
low temperatures (f7|). 
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